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A  CRITICAL  EVALUATION  OF  NEW  PLATE  THEORIES 


APPLIED  TO  LAMINATED  COMPOSITES* 


Charles  W.  Bert 

School  of  Aerospace,  Mechanical  and  Nuclear  Engineering 
University  of  Oklahoma 
Norman,  Oklahoma 


ABSTRACT 


The  plate  theory  recently  developed  by  Levinson  is  extended  to 
laminates.  Closed- form  solutions  of  this  theory ,  as  'Sell  as  those 
of  Reissner-’lindlin  plate  theory  with  appropriate  shear  correction , 
Seide’s  discrete-layer  plate  theory ,  and.  Lo,  Christensen,  and  Wu's 
higher-order  theory  are  all  compared  with  Pagano’s  elasticity- 

C 

theory  solution  for  the  cases  of  cylindrical  bending  of  a  single 
orthotropic  layer  and  a  symmetric  cross-ply  (C°/00°/0° )  laminate 
consisting  of  three  ecp.ial-thickness  layers.  Quantities  compared 


are  maximum  plate  deflection,  bending  stress 
transverse  shear  stress  distribution. 


distribution,  and 


Dist 


1  INTRODUCTION 

It  has  long  been  known,  through  Saint-Venant' s  flexure  theory  as  well  as 
through  experimental  observations,  that  the  elementary  Bernoul 1 i -Euler  beam 
theory  is  inaccurate  except  in  the  case  of  pure  bending  (no  transverse  shear 
forces)  or  very  slender  geometry  (large  length/depth).  Mindlin  and 


*  An  abbreviated  and  preliminary  version  of  this  paper  was  presented  at  the 
Symposium  on  Mechanics  of  Composite  Materials,  sponsored  by  the  Applied 
Mechanics  Division,  American  Society  of  Mechanical  Engineers  Winter  Annual 
Meeting,  Boston,  MA,  Nov.  13-18,  1983. 
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Deresiewicz1  provided  an  excellent  review  of  early  improved  theories  to  take 
into  account  transverse  shear  deformation,  including  the  pioneering  work  of 
Bresse^  in  1859  and  analogous  work  by  Timoshenko  in  1921-22. 3,4 

A  similar  situation  was  evident  in  the  theory  of  plates,  since  the 
classical  thin  plate  theory  (CPT)  due  to  Germaine  and  Lagrange  suffers  from 
the  same  deficiencies  of  Bernoulli-Euler  beam  theory,  namely 

a)  Transverse  shear  strain  is  neglected 

b)  In-plane  normal  strain  is  distributed  linearly  through  the  thick¬ 
ness,  rather  than  nonlinearly 

c)  Transverse  normal  strain  is  neglected 

Reissner5,e  and  Mindlin  presented  generalizations  of  the  Bresse- 
Timoshenko  beam  theory  to  plates  and  thus  made  the  first  attempts  to  include 
transverse  shear  deformation  in  plate  theory.  Their  theories  differed  not 
only  in  application,  References  5,6  to  static  problems  and  Reference  7  to 
dynamic,  but  also  in  the  definition  of  the  kinematic  parameters.  These 
theories  not  only  suffered  from  deficiencies  b  and  c,  listed  above,  but  also 
the  transverse  shear  strain  was  distributed  linearly  through  the  thickness. 
This  required  use  of  a  transverse  shear  correction  factor,  either  implicitly 
(Reissner)  or  explicitly  (Mindlin). 

Perhaps  the  first  attempts  to  account  for  a  more  realistic  distribution 
of  transverse  shear  strain  were  due  to  Ambartsumyan  in  1957;  see  Reference  8, 
page  40.  In  his  theory,  all  deficiencies  (a,b,c)  were  removed.  Later  work 
conducted  in  this  same  spirit  was  due  to  Reissner9  and  Levinson10,  and  for 
the  geometrically  nonlinear  case,  Schmidt11.  Numerous  so-called  high-order 
theories  were  presented  by  Donnell12,  Tiffen13,  Tiffen  and  Sayerlu,  Tiffen 
and  Lowe15,  Lee16,  Berdichevskii17,  Pane13,  Lo  et  al.19,  Cheng20,  Celep21, 


Krenk22,  Voyiadjis  and  Baluch23,  and  Shirakawa24.  An  interesting  tabular 
comparison  of  various  plate  theories  was  presented  by  Irretier". 

The  first  complete  laminated  anisotropic  plate  theory  is  generally 
attributed  to  Reissner  and  Stavsky2f;  this  was  a  laminated  version  of  CPT. 
However,  it  has  been  known  for  a  long  time27  that  transverse  shear  deforma¬ 
tion  plays  a  considerably  larger  role  in  structures  made  of  filamentary 
composite  materials  than  those  of  isotropic  materials.  The  explanation 
for  this  is  the  very  low  shear  moduli,  relative  to  in-plane  elastic  moduli, 
exhibited  by  such  composite  materials.  The  laminated  plate  versions  of 
the  Mindlin  and  Reissner  plate  theories  are  due  to  Yang  et  al.2'-  and  Whitney 
and  Pagano:°,  respectively.  In  Reference  28,  a  dynamic  method  of  determin¬ 
ing  the  appropriate  shear  correction  factors  was  introduced  and  in  References 
30  and  31,  static  methods  were  introduced.  Also,  the  works  of  Reissner''  ’ 3  ‘ 
Bondar' and  Rasskazov34,  and  Green  and  Naghdi  r  should  be  mentioned. 

Higher-order  laminated  plate  theories  were  reviewed  by  Lo  et  al.3f', 
who  applied  the  laminated  version  of  their  own  high-order  theory1'.  Also, 
the  work  of  Whitney  and  Sun37,  Librescu33,  Murthy5\  Soni  and  Pagano°,  and 
Rehfield  and  Vallisetty1’1  should  be  mentioned. 

It  is  noted  that  in  all  of  the  laminated  theories  discussed  above,  the 
shear  angle  was  either  not  permitted  to  vary  at  all  from  layer  to  layer  (in 
the  theories  of  References  28  and  29)  or  required  to  vary  in  a  smooth,  a  priori 
fashion  in  the  higher-order  theories.  Apparently,  the  first  attempts  to  con¬ 
sider  each  layer  in  a  laminate  as  a  separate  beam  or  plate  are  due  to  Refer¬ 
ence  42  for  the  static  case  and  Reference  43  for  the  dynamic  case  of  multi¬ 
core  sandwich  beams+and  Reference  44  for  laminated  beams  loaded  statically. 

+A  sandwich  beam  is  generally  understood  to  be  one  having  two  or  more  relatively 
stiff,  thin  layers  (called  facings)  and  one  or  more  relatively  flexible,  thick 
layers  (called  cores) . 
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This  approach  was  also  used  by  Ambartsumyan8,  page  75,  for  thick  plates, 
considering  each  layer  as  an  Ambartsumyan-theory  plate.  See  also  the  shell 
work  of  Hsu  and  Wang45  and  the  plate  work  of  Seide46,  who  considered  each 
layer  as  a  Reissner  plate. 

Here,  the  terminology  "smeared  laminate  model"  (SLM)  is  used  to  des¬ 
cribe  laminate  theories  of  the  type  (References  28-29,  32-36);  "discrete 
layer  model"  (DLM)  is  used  to  describe  laminate  theories  of  the  type 
(References  42-46). 

2  METHODOLOGY  OF  COMPARISON 

The  objective  of  an  improved  theory  for  laminated  plates  is  to  achieve 
greater  accuracy  of  prediction  than  is  possible  with  classical  thin  plate 
theory  (CPT)  or  even  classical  (Reissner-Mindl in)  shear  deformable  plate 
theory  (SDT),  without  requiring  the  complexity  of  three-dimensional  elasti¬ 
city  theory  or  even  that  of  the  more  complicated  higher-order  plate  theories. 
It  is  customary  to  evaluate  the  accuracy  of  various  improved  theories  by 
comparison  of  the  results  for  a  specific  situation  with  those  of  a  three- 
dimensional  elasticity  solution  for  the  same  situation.  In  this  regard, 
the  closed-form  solution  due  to  Pagano47,  for  cylindrical  bending  of  a 
simply  supported  laminate  under  a  sinusoidally  distributed  normal  pressure 
loading,  has  been  widely  used. 

In  the  present  work,  four  different  theories,  described  in  the  ensuing, 
are  applied  to  cylindrical  bending  in  two  different  cases: 

Case  1  Homogeneous  orthotropic  material 

Case  2  Symmetric  cross-ply  laminate  (three  layers) 

In  both  cases,  the  material  considered  is  the  same  as  that  of  Pagano1*7: 
E^/Ey  =  25,  E|_/G[^-|-  =  50,  ^y/Gyy  =  2.5,  V|^y  =  Vyy  =  0.25. 


This  theory  is  due  to  Yang,  Norris,  and  Stavsky:8  and  to  Whitney  and  Pagano-9. 
It  is  used  here  in  conjunction  with  the  general  shear  correction  factor  de¬ 
rivation  given  in  Reference  31  or  more  explicitly  in  Reference  48. 


Classical  Shear  Deformable  Plate  Theory  (Discrete  Layer  Model) 

This  theory  is  due  to  Seide46,  who  applied  it  and  worked  it  out  in  detail 
for  a  symmetric  three-ply  laminate. 

Laminated  Version  of  Levinson's  Theory  (Smeared  Laminate  Model) 

This  theory  is  presented  in  Appendix  A.  Due  to  the  nonlinearity  of  the 
axial -normal -stress  distribution  through  the  thickness,  use  of  the  equilibrium 
equation  for  the  xz  plane  requires  a  higher  degree  of  nonlinearity  in  the 
shear-strain  distribution  than  is  assumed  a  priori  in  this  theory.  Thus, 
contrary  to  the  remarks  in  Reference  19,  a  shear  correction  factor  may  still 
be  needed  in  this  theory.  However,  due  to  the  complexity  of  the  shear-strain 
distribution  resulting  and  the  dependency  of  k  upon  the  normal  pressure  (see 
Appendix  B),  it  is  not  practicable  to  work  out  this  correction  factor  in  general 
However,  it  is  shown  that  for  the  homogeneous  case  and  p  =  0  that  k  =  1. 

Lo-Christensen-Wu  (LCW)  Higher-Order  Theory  (Smeared  Laminated  Model) 

This  theory  was  presented  in  Reference  19  for  the  homogeneous  case  and  in 
Reference  36  for  the  laminated  one. 


4  RESULTS  AND  DISCUSSION 

Homogeneous  Case 

In  this  case,  due  to  the  absence  of  bending-stretching  coupling,  B  =  C  = 


and  furthermore,  E/D  =  1/5.  Then,  eqn  (A-15)  gives  for  the  Levinson  theory 


“nax  5  OV0-"*  *  <6/S)<p0/QS5h»  )  (1) 

This  result  is  identical  to  that  of  classical  Reissner  type  theory'’,  i.e., 
Whitney-Pagano29  with  a  shear  correction  factor  k  =  5/6. 

For  comparison  with  Pagano's  numerical  results1*7,  eqn  (1)  can  be  re¬ 
written  as  follows: 


w  =  100  ETh’wmax/Poeu  =  (100/-‘*)[12(1 -vLTvTL)  +  (6~75)(EL/GLT)(h/0  ](Et/El) 

(2) 

Using  the  previously  mentioned  material -property  ratios  in  eqn  (2),  one 
obtains  w  =  1.981  for  h/2  =  1/4,  considered  by  the  present  investigator  to  be 
the  maximum  thickness  of  a  plate  rather  than  a  block.  This  value  compares 
very  favorably  with  a  value  of  approximately  1.95  read  from  the  curve  of 
Reference  47.  To  the  small  scale  of  the  plot  in  Reference  19,  this  is  in 


agreement  with  the  LCW  higher-order  theory. 

Again  for  comparison  with  Reference  47,  eqn  (A- 1 6 )  can  be  used  to  obtain 


<3x>max  5  <°*Wo  ■  («/’->(*/*•>  *  (1/10){El/Glt)(1  -  yj_jvTL) 


For  h/£  =  1/4,  eqn  (3)  gives  (o  )  =  14.74,  which  is  in  fairly  good 

x  ma  x 

agreement  with  Pagano's  elasticity-theory  value  (approximately  14.1).  The 
prediction  of  classical  Reissner  SDT  is 

<5x>max  ’  (W-)U/h)2  (4) 

which  yields  a  value  of  (ov)  =  9.727,  which  obviously  is  considerably 

x  max 

inaccurate.  Incidentally,  Fig.  2  of  Reference  19  appears  to  be  drawn  in¬ 
accurately  in  the  vicinity  of  h/2  =  1/4  and  thus  cannot  be  used  to  compare 
the  LCW  value  with  the  above. 


Three-Layer  Cross-Ply  Laminate 


Again,  the  absence  of  bending-stretching  coupling  causes  coupling  stiff¬ 
nesses  B  and  C  to  vanish.  Also,  it  can  be  shown  that 


D  =  (26  El  +  ET)(h3/324)(l  -  vltvtl)-1 

E  =  (242  El  +  Ey)  (h  Vl4,580)  (1  -  v^)'1,  S55  =  (23GLT  +  26GJT)  (h/81 ) 


Then,  application  of  eqn  (A-15)  gives,  for  h/f.  =  1/4, 

*  100  EThWv4  =  2-630 

This  deflection  value  is  approximately  12.3%  lower  than  the  exact  value  of 
approximately  3.0  obtained  in  Reference  47. 

To  apply  SOT,  the  following  equation  for  the  shear  correction  factor 
is  used  (References  31,48) 


k  =  D/I 


f  G(k)dz 


[b'/G^]  dz) 


(6 


or,  for  the  specific  laminate  (three  equal -thickness  plies) 

l  _  5  _  (26e  +  1  )2g  _ 

'  6  (1  +  2g ) [ 1 02e-  +  g(120e~  +  20e  +  TJ1 


(Reference  48) 


(7; 


where  e  s  E^/Ey  and  g  £  G^y/Gyy.  For  e  =  25  and  g  =  2.5,  eqn  (7)  yields 
k  =  0.5828.  Also,  and  D  are  given  by 


G(k)  dz  =  (2  Ga  +  GTT) (h/3)  ,  D  =  J  z  Qn(k)dz 


26El+Et  ^ 
324(1  -  vLTvTL) 


Then,  the  maximum  deflection  is  given  by 


(8) 


Wmav  “  (Pn/8i4)  +  (p  /kSf-ca") 


(9 


The  result,  for  h/t  =  1/4,  w  =  3.226,  which  is  7.5%  higher  than  the  exact  one. 
A  comparison  of  the  present  deflection  results,  for  h/i.  =  1/5,  with  those  of  a 
number  of  other  investigators  is  given  in  Table  1.  The  bending-stress  and 
shear-stress  distributions  are  shown  in  Figures  1  and  2. 

5  CONCLUDING  REMARKS 

On  the  basis  of  the  comparisons  made  in  this  study,  it  is  concluded  that  the 
Levinson-type  theory  is  more  accurate,  i.e.,  closer  to  the  exact  elasticity 
solutions,  than  classical  SDT.  The  higher-order  LCW  theory  is  also  more 
accurate,  but  requires  too  much  computation  to  justify  the  accuracy  achieved. 
Both  of  these  theories  predict  the  nonlinear  distribution  of  bending  stress 
through  the  thickness,  while  SDT  as  well  as  CPT  does  not. 

However,  the  Seide  theory,  which  is  a  discrete  layer  version  of  classical 
SDT,  is  more  accurate  in  predicting  shear-stress  distribution  than  any  of  the 
smeared  laminate  theories  mentioned  above.  Unfortunately,  the  Seide  theory 
is  not  quite  as  accurate  in  predicting  the  maximum  bending  stress,  since  it 
predicts  a  linear  distribution  of  bending  stress.  This  suggests  that  a  new 
theory,  a  discrete  layer  version  of  the  Levinson  theory,  should  be  most 
accurate.  It  would  be  expected  to  have  an  accurate  prediction  of  shear-stress 
distribution  like  the  Seide  theory  and  an  accurate  prediction  of  the  nonlinear 
bending  stress  distribution  like  the  Levinson  and  LCW  theories. 
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APPENDIX  A:  LAMINATED  VERSION  OF  LEVINSON  PLATE  THEORY 


Taking  Cartesian  coordinates  x  and  y  in  the  plane  of  the  plate  and  z  as 
the  thickness  normal  coordinate,  measured  positive  downward  from  the  midplane 
of  the  laminate,  one  starts  with  the  following  displacement  field 


U(x,y,z)  =  u  (x,y)  +  zc  (x,y)  +  z3d  (x,y) 

o  x  x  (A.-, 

V(x,y,z)  =  vQ(x,y)  +  zoy(x,y)  +  z?$y(x,y)  ;  W(x,y,z)  =  w(x,y) 

It  is  noted  that  the  midplane  displacements,  not  included  in  Levinson's 
original  theory10,  are  necessary  here  in  order  to  provide  for  the  bending¬ 
stretching  coupling  exhibited  in  unsymmetric  laminates.  Further,  as  was 
pointed  out  by  Murthy?~‘,  terms  in  z:  are  not  needed  due  to  the  requirement 
of  zero  shear  stress  (and  thus  zero  shear  strain)  on  the  upper  and  lower 
surfaces  of  the  laminate. 

The  thickness  shear  strains  are 


'xz 


(A-2) 


where  (  )  =  d{  )/3x,  etc. 

Imposing  the  condition  of  zero  shear  strain  at  the  laminate  surfaces 
(z  =  +  h/2),  one  can  express  t  and  $■  in  terms  of  some  of  the  other  kinematic 

X  y 

variables; 

*x  =  - (4/3h2 ) (c  +  w  )  ,  i  =  - (4/ 3h~ ) (y  +  w  J  (A-3) 

A  A  j  J  yj 

Thus,  the  strain  field  is 


U  +  Z'i 
o,x  vx,x 


=  V 


V  +  z- 

n  yj 


(4/3h  2)23(<^XjX  +  w^xx) 
(4/3h2  )z3(i''  w  +  w  ) 


where 


544  S45 

545  S55 


:v  +  w  vl  ( i . j=l  ,2,6) 

j 


+  w 
X  ,x 


h/2 

(Aij-Bij,Dij>  '  J  dz 

-h/2 


(Ci  j  * E -j  j )  =  C4/3h  ) 


h/2 


h/2 


■h/2 


(1,2)2  Q..  dz 

*  J 


skt=|  [1  -  (2z/h)']Qkl  d2 

-h/2 


>  (l.j-1 .2,6) 


k ,  =4 , 5 


The  equilibrium  equations  of  elasticity  are 


°i j »j  +  Fi  ‘  0 


Integration  of  eqns  (A-9)  through  the  thickness  of  the  laminate  and 
of  eqns  (A-6)  yield  the  usual  plate  equilibrium  equations 

N  +  N  +  p  =  0 
x,x  xy,y  Kx 

N  +  N  +  p  =  0 
xy,x  y,y  Hy 

V  +  V  +  p  =  0 
x.x  y,y  w 


1  +  M 

V 

+  m 

=  0 

x , x  xy,y 

X 

X 

1  +  M 

V 

+  m 

=  0 

xy,x  y,y 

y 

y 

where  p.  and  m^  are  body  forces  and  body  moments,  and  p  is  the  normal 


(A-7) 


pressure. 


Substitution  of  the  laminate  constitutive  eqns  (A-7)  into  eqns  (A-10) 
yields  finally  five  plate  equilibrium  equations  in  terms  of  the  five  gener¬ 
alized  displacements  u  .  v  ,  w  .  v  .  C-  . 

o  o  o  x  y 

For  the  case  of  cylindrical  bending,  all  derivatives  with  respect  to  y 
vanish  and  the  equilibrium  equations  become 


A, ,u  +  B(,;'  -  CnW 

11  o,xx  11  x,xx  11  ,xxx 

S55^x,x  +  w,xx ^  P  =  0 

Blluo,xx  +  Dllvx,xx  '  Ellw,xxx 


=  0 


-  S(;  +  w  )  =  0 

X  ,x 


( A- 1 1 ) 


These  can  be  uncoupled  to  yield 


(D11  ‘  A^  w , xxx  =  p  '  (D11 


Jl 

Si 


Bncn 

An 


E11)p,xx/S55 


( A- 12) 


For  a  sinusoidally  distributed  pressure 
p  =  pQ  sin  ,ix  ;  n  =  -/z 

the  deflection,  for  freely  supported  edges,  is 

w  =  wmav  sin  ax 
max 


where 

_  .  po _  ,  _po  ,,  .  BC  -  AEX 

Wmax  "  [D  -  (B  /A)]a4  U  AD  -  B‘ ' 


( A- 1 3) 


(A- 14) 

(A- 15) 


and  the  subscripts  11  have  been  omitted  from  A,  B,  C,  D,  and  E  for  brevity. 

In  eqn  (A-15),  the  first  term  is  identical  to  that  of  CPT  and  the  quantity 


1  +  (BC  -  AE) (AD  -  B2) 

is  the  multiplier  of  classical  shear  deformable  plate  theory. 

It  is  interesting  to  note  that  the  normal  strain  is  not  distributed 
linearly  through  the  thickness  but  has  a  term  in  the  cube  of  z: 

£,  -  [(B/A)  -  z]w  -  fc!  +2-4  rr)(P/Scc) 


(A-16) 


APPENDIX  B:  CORRECTION  COEFFICIENT  FOR  A  LAMINATE  IN  LEVINSON  THEORY 
There  are  numerous  methods  of  determining  the  shear  correction  coefficient, 
some  dynamic  and  some  static.  The  two  most  popular  static  methods  are: 

1.  Through  use  of  the  axial -force-equil ibrium  equation  of  elasticity. 

This  approach  is  identical  to  that  of  Jourawsky's  theory  of  transverse  shear 
in  beams  (and  Reissner's  for  plates)  and  yields  a  value  of  5/6  for  a  homo¬ 
geneous  rectangular  section. 

2.  Through  use  of  equivalency  with  the  results  of  Saint-Venant  theory 
of  flexure.  This  method,  originated  by  Cowperi,:‘,  yields  a  value  dependent 

upon  Poisson's  ratio  for  a  homogeneous  rectangular  section  of  isotropic  material 
Here,  the  former  approach  is  used.  The  plane-strain  equilibrium  of 
axial  forces  in  the  xz  plane  is  expressed  by 


o  +i  =  0 

x,x  xz,z 


Thus, 


xz 


f 

! 

-h/2 


a  v  dz 
x  ,x 


(B-l ) 

(B-2) 


However, 


=  x  ,/Q 


xz 


(k) 
xz,w55 


(B-3) 


.th 


where  the  superscript  k  denotes  the  k  layer. 
Thus, 


v  ■ 


-h/2 


Q<k)  E 
11  x,x 


dz 


(B-4 ) 


or 


’xz  ■ 


■h/2 


Q,(k)[u  +  zv  -  (4/3h2 )z 3 (;.  +  w  )]dz  (B-5) 

yll  L  o,xx  vx,xx  '  x,xx  ,xxx/J  ' 


Now,  it  is  necessary  to  express  the  generalized  displacement  in  terms  of 


the  generalized  forces  by  inversion  of  the  laminate  constitutive  relations, 
eqns  (A-7).  The  results  are 


DN  -  BM  nr 

u  =  _ *  +  2L 

uo,x  AD  -  B-  AD 


AMx  -  BNx  ,  AE 
Vx  "  AD  -  B2  AD 

i i>  +  w  =  V  /Scc 
x,x  ,xx  x  55 


“F  (Vx,x/S55^ 


^F  {Vx,x/S55> 


(B-6) 


where  subscripts  11  have  been  omitted  from  A,  B,  C,  and  D  for  brevity. 
Substituting  relations  (B-6)  into  eqn  (B-5),  one  obtains 


n.0(k),  f2  0(k),0Nx.x-BM*.x*  (CD-BEHVx>xx/S55) 
\z  -  '[1/Q55  ]  Q11  1 -  AD  -  B 


-  bnx,xMae-bc)(vx-xX/s55) 


AO  -  B- 


Now,  beam-type  equilibrium  requires  that 


z  -  (4/3h')(Vx>x/S55)z  }dz 


(B-7 ) 


N  =0  ,  M  =  V  ,  V  =  -  p 

X,X  X,X  X  X ,x 


(B-8) 


Thus,  eqn  (B-7)  reduces  to 


■  -tVog’nfcft  Vx  +  <CP^S55>  +  [L'D:'BEi-AD*-(r:-B-Clb-](-P,x/S55)I 


(B-9) 


where  the  partial  stiffnesses  are  given  by 


(a,b)  z  j  (l,z)Qj*;  dz 
-h/2 

z 

c  =  (4/3h2)  z3q{^  dz 


(B-10) 


For  a  symmetric  laminate,  B  =  C  =  0  and  eqn  ( B- 1 0 )  reduces  to 
■'xz  “  -  '/Q^>C(b/D)Vx  ♦  (C/S55)p  -  (Eb/DS55)p>x] 

To  derive  the  shear  correction  factor  (k),  one  uses 


(B-l  1 


k  U$C  •  U$E 


(B-l  2 


C  E 

where  and  are  the  shear  strain  energies  (per  unit  length)  calculated 
on  the  basis  of  the  constitutive  relation  and  the  equilibrium  equation  of 


elasticity,  respectively.  Thus, 


h/2  :  h/2 

usc  =  H  o  -  *  s  dz 


(B-l  3 


where 


+  W,x  *  -<0  •  :x  =  -(4/3h'  >  ’o 


(B-14 


A1  so, 

h/2 

USE  4  I  1  t(W°)Vx  +  -  (eb/DS55)P  xr  dz 

-h/2 


Vx  =  S55vo 


Thus,  eqn  (B-15)  can  be  rewritten  as 


USE  *  ?  j  ^[(bS55/D)^o  *  (C/S55)p  -  (Eb/DS55)p>x]£<lz 


Finally,  in  view  of  eqn  ( B-l 2 ) 


(B-15 


(B-16 


( B-l  7 
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h/2 

-h/2 


(1  -  4  p)  G  dz 


h/2  2 

G  ^bS55^Yo  +  ^C//S55^p  "  ^Eb/DS55^p,x^  dz 

-h/2 


(B-18) 


It  is  noted  that,  in  contrast  to  the  shear  correction  factor  in  Bresse- 

Timoshenko-Reissner  theory  which  is  independent  of  p,  the  present  expression 

for  k  depends  upon  p  and  p  .  Since  for  the  sinusoidal  distribution  of  p, 

*  x 

the  quantities  p  and  p  vary  with  x,  eqn  (B-18)  implies  that  k  must  vary 

i  x 

with  x.  However,  it  is  not  practicable  to  consider  this  variation  here. 

Thus,  for  p  =  p  =0,  eqn  (B-18)  reduces  to 
*  X 


h/2 

(1  -  4  p)  G  dz 

_ 

k~  *  -  (B-l 9) 

h/2 

(S55/0)'’  |  (b:/G)  dz 

-h/2 


For  the  homogeneous  (single-layer)  case,  and  G  are  independent  of 

z  and 

b  =  (l/2)[z2  -  ( h '  / 4 ) ]  Qn 


Then,  eqn  ( B-l 9)  gives  k  =  1  precisely. 


TABLE  1 


Dimensionless  Deflection  of  a  0°/90°/0°  Laminate  by 


Various  Theories,  for  h/e  =  1/5 


Theory 

Reference 

w  5  100ETh7w  /pc.'* 

T  max  ^o 

Exact 

47 

%  2.0* 

Seide 

46 

%  2.0* 

Murthy 

39 

%  1.87* 

Present  (SDT) 

• 

2.25 

Present  (Levinson) 


1.87 


DIMENSIONLESS  BENDING  STRESS,  a  i,z)/p 


Fig.  1  Bending  stress  distribution  through  half  thickness  of 
three-ply  laminate;  h/£  =  l/4.  Stress  in  middle  ply 
too  small  to  show  at  this  scale. 
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